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1 Introduction 



There are very few exactly solvable spin systems, in the sense that their complete spectra can 
be determined in closed form. The most notable one is the Ising model in transverse field or 
the XY model [1], where the spectrum can only be determined with the help of highly nonlocal 
Jordan- Wigner transformation, rendering the calculation of correlation functions very difficult. 
For the other example, the Heisenberg spin chain, the energy eigenvalues and eigenvectors can 
only be determined implicitly, by solving the system of coupled non-linear Bethe ansatz equa- 
tions [2], [3] . For a large class of spin models, only the ground state can be found by the matrix 
product approach [H El [6l [H El [3 |10] . Having an exactly solved quantum many body system, 
in the sense of complete determination of its spectrum, is always a fortunate situation, which 
enables one to make a detailed study of thermal and dynamical properties. It also allows one 
to solve other related systems by perturbation techniques. 

While the study of many body systems has been traditionally done in the community of 
condensed matter and statistical physics, and also by mathematical physicists interested in ex- 
act solutions, in recent years, these systems have attracted a lot of attention from the quantum 
information community. The reason is at least two-fold: on the one hand an array or lattice 
of two-level quantum systems (qubits) is the natural candidate for implementation of quantum 
information processing tasks and on the other hand, concepts and tools developed in quantum 
information [11], which are mainly aimed at characterizing the nature of quantum states like 
their entanglement, have been quite useful in understanding different phenomena in such sys- 
tems, i.e. quantum phase transitions [T2| [13]. 

While in condensed matter physics, the focus is on the Hamiltonian and interactions, in 
quantum information, the emphasis is on the quantum states and their properties, i.e. their 
bi-partite and multi-partite entanglement. Needless to say, we are not always in the happy situ- 
ation to have both a physically plausible Hamiltonian on the one hand and an easily obtainable 
spectrum on the other. For example in the matrix product formalism, although we can con- 
struct ground states with desired symmetries, it is not guaranteed that the parent Hamiltonians 
are of real experimental interest. Fortunately with the recent advances in optical lattices and 
cold atoms, we have more freedom in manipulating systems of many body two-level systems 
(qubits). Therefore there is less reservation than before in proposing many-body Hamiltonians 
whose ground states, or the low level excited states, may have desirable properties for imple- 
mentation of quantum information processing. 

In this regard, an interesting model has been recently proposed [H] which has the very 
desirable feature of showing topological order and anyonic excitations. This was the first model 
for topological implementation of quantum computation, where qubits are encoded into the ho- 
mological classes of loops on a surface which are hence resistent to local errors. The important 
point is that topological order is not related to the symmetries of the Hamiltonian and it is 
robust against arbitrary local perturbations, even those that destroy all the symmetries of the 



1 



Hamiltonian. 

On any lattice, the Kitaev model is defined by the Hamiltonian 

H=-jY.As-KY.Bp, (1) 

s p 

where J and K are positive coupling constants, and the vertex operators Ag, and plaquette 
operators Bp are constructed from Pauli operators X = ax and Z = az a.s follows: ^4^ is the 
product of all X's on the links shared by the vertex s and Bp is the product of all Z operators 
around a plaquette. These operators commute with each other for any geometry of the network. 

It is well-known that the ground state of this model has symmetries not inherited from the 
Hamiltonian but from the topology of the surface on which the model is defined. On a genus 
g surface without boundary, the ground state of model ([I]) has 4^ fold degeneracy with a gap 
which cannot be removed by local perturbations. Therefore on such a surface, the ground 
space can encode 2g qubits, in a way where resistance to errors is automatically ensured by 
the topology. In recent years, there has been intensive activity on this model and its variations 
and generalizations in many directions, see the works [IHl UHl [171 HH |19l EQl EU [221 ESI l2l| and 
references therein for a sample. The ground states of the Kitaev model can be characterized 
rather simply in a formal way. Since = B^ = I, the ground state is one which is stabi- 
lized by all these operators, i.e. As\g.s) = Bp\g.s) = \g.s). One such state can be written as 
\g-s) = J2cCz\^), where \Q) is a sea of spins in +x direction and Cz ■= Ylt^c is a product 
of flipping operators Z around the closed loop C. If C is a homologically trivial loop, i.e. if it 
is the boundary of a region on the surface, then Cz can be expressed as a product of Bp oper- 
ators. Since in \g.s) we are summing over all such loops, the action of any Bp on \g.s) leaves it 
invariant. Since any As commutes with all Bp operators and A^lfi) = 1^2), it is immediate that 
As\g.s) = \g.s) and hence \g.s) is actually a ground state of ([1]). The other degenerate ground 
states are obtained from \g.s) by acting on it by product of flipping operators Zi around non- 
trivial homology cycles (Closed loops which are not boundaries of regions are called homology 
cycles or just cycles. On genus one surfaces these are the same as non-contractible loops, but 
on higher genus surfaces they are different.) 

The excited states are formed by enacting on the ground states by flipping operators along 
open strings and hence creating two particles, [25] called anyons due to their exchange prop- 
erties. While the ground state of ([1]) is rather simple and in fact its entanglement properties 
and that of the related topological color codes [26] have been studied in a number of works 
[271 [28], a complete characterization of the spectrum is difficult, due to the exponentially large 
number of open string configurations. In other words, it is known that any collection of open 
strings create an excited state, but determining the degeneracy of such excited states is not 
simple, due to the above-mentioned difficulty. 

One may expect that on a lattice with a simpler structure, these problems can be overcome. 
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In this regard, spin ladders may be of interest, not only due to their own interest as systems 
interpolating between one and two dimensional systems, but also since they can be used for 
approximate solution of the more physical two dimensional systems, when the latter can be 
approximated as an array of ladders with negligible couplings between them. 

In fact the simple structure of ladders facilitates the study of many interesting phenomena 
which are otherwise difficult to study in general lattices. For example one of these phenomena 
is the dynamics of defect production, in passing a critical point when a phase transition occurs 
[29] . This is usually captured in what is known as the Kibble-Zurek scaling law. However it has 
recently been shown [30], through a detailed study of a ladder system, i.e. the Creutz ladder 
[31] , that in systems with topological order, edge states can dramatically modify this scaling law. 

It is the purpose of this paper to make a complete study of the Kitaev model on spin ladders. 
The ladder with periodic boundary condition in one direction has the topology of a cylinder and 
its ground state is doubly degenerate. We will determine the energy spectrum completely and 
from there we calculate the partition function and the total entropy. Then we will determine 
the reduced density matrices and entropy of various subsystems, where it is found how the 
cylinder topology affects these properties. Finally we discuss briefly how to extend the method 
to three and higher leg ladders and eventually to the two dimensional square lattice. 

The structure of this paper is as follows. In section ([2]) we explain the model and in its 
various subsections, calculate the spectrum, the partition function and the thermal averages 
of non-local operators. In section ([3]) we calculate the density matrices and entropy of various 
subsystems which show among other things that there is no entanglement between any two 
spins at any temperature. We end with a discussion on how to extend the method to larger 
lattices. 



2 The Kitaev Model on the ladder network 

Consider a two-leg ladder network with length with the labeling of links as shown in figure 
([T]), where we use the superscripts -|- and — to denote the vertices pertaining to the lower and 
upper legs. For this ladder, the vertex and plaquette operators take the following form, 

Af := Xi/_iXiXii , A^ := Xi"^iXiXi" , Bi = ZiZi^iZi^Zi" , (2) 

Throughout the paper, we use the generic names Ag and Bp for pointing to general vertex and 
plaquette operators in a network and the names indicated in Eq. ([2]) for pointing to specific 
operators in the ladder. For a closed surface the operators in ([1]) are constrained by the relations 
Os As = Tip Bp = I, however for the ladder, which is a surface with two boundaries, the second 
constraint does not hold and we are left with 

IlAs = L (3) 
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Figure 1: The labeling that will be used in the text for the links on the ladder. 

Therefore the operators in ([2]) are 3A^ — 1 commuting operators in the 2'^^ dimensional Hilbert 
space of the ladder. They all commute with the Hamiltonian and hence the ground state is 
two-fold degenerate. The following non-local operators play an important role. 

W,:=\{Z, l^, :=n^n (4) 

i' i 

It can be clearly seen that they commute with each other and also with the Hamiltonian, 

[W,, W,] = [H, W,] = [H, W,] = 0. (5) 

They correspond to a cycle around the ladder. Both of them square to 1 and hence have 
eigenvalues ±1. As we will see the operator Wz generates the two-fold degeneracy of the 
ground state. 

2.1 The ground and the top states 

It is clear from ([T]) that the ground state is a state with the property 

A.l^o) = l^o), B,\^o) = l^o), Vs, I. (6) 

To find the explicit form of the ground state, we use the property Bf = I which is equivalent 
to Bi{l + Bi) = 1 + Bi and introduce the state 

\^o)■■=^m^ + B^)\n+), (7) 

where |f2+) := 1+)*^^^, is a product of all spins in the |+) := \x,+) direction. Note that 
X|±) = ±|±) and Z\±) = |=i=). It is obvious that l^'o) satisfies the condition ([6]) and hence 
is a ground state. The other ground state is obtained from |\l/o) by the action of the operator 
Wz which corresponds to the only nontrivial homology cycle of the surface, that is |\l/o) and 
l^o) '■— W^l^o) form the doubly degenerate ground states of the model, 

Hi\^,), \%)) = -N{2J + mm,), \%)). (8) 

Note that Wz, being a cycle, cannot be expressed as any combination of product of plaquette 
operators Sj, which are all homologically trivial. This shows that the two states |\l/o) and |\1/q) 
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are independent. 



The top state, the state with the highest energy is one which is an eigenstate of all the vertex 
and plaquette operators with eigenvalues —1. It can be readily verified that the following state 
is such an state: 

\'^top):=^W-B,)\n^), (9) 

where = |— ^i^^Af jg ^^iq product of all spins in the |— ) := \x, — ) direction in the network. 
The other degenerate state is obtained by the action of the operator W^, i.e. 1^^^^) = Wzl'^top)- 
The highest energy will be = N{2J + K), 

Hilltop), l^'top)) = N{2J + K){\m,,,), \%^;)). (10) 



2.2 The complete spectrum 

We can now construct the full spectrum. To this end we note that for any arbitrary state 
the state, 

\{{l + {-lfB,)\x), k= 0, 1, 

i 

is an eigenstate of all the plaquette operators, Bj with eigenvalues (—1)'^. The reason is the 
relation Bj{l ± Bj) = ±{1 ± Bj). Hereafter we use the abbreviation Bi for such a string of 
operators, 

■■= 11(1 + i-iy^B,). (11) 

i 

These operators satisfy 

BiBi, = 2^<5i,r5i. (12) 

To construct states which are eigenstates of the vertex operators and at the same time be 
independent, let us define the following operators 

Ar,s:=l[ZrZp, (13) 

i 

where the sequence of and Sj are or 1 and the labeling are those shown in figure ([T]). Figure 
([6] ) shows the links which contribute to the construction of such operators. It is important to 
note that the links of only one leg are among this set. Now let us define 

l^i.r.s) := Ar,s|^i) := -^ArMn+). (14) 



The operators in f|T3l) have simple commutation relations with the vertex and plaquette oper- 
ators. One can verify the following relations 

BjAr^s = Ar^sBj 
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A+Ar,s = (-l)^^+^-^+^^Ar,sA+ (15) 

This leads to 

Bj\^l,r,s) = (-l)'^|^l,r,s) 
^7l^I,r,s) = (-l)^^|^l,r,s) 

^;i^I,r,s) = (-l)^^+^-^+^^|^l,r,s). (16) 

The above relations indicate that the states thus constructed are eigenstates of Hamiltonian, 

^|^l,r,s = £^l,r,s|^l,r,s) (17) 

where the energy is found from ([T]) and fll6p to be 

^l,r,s(^, i^) = - E [Jii-^Y' + (-l)'-^+^-^+^0 + A'(-l)'^] . (18) 

j 

In terms of Ising-like variables, Ri := {—1Y\ Si := (—1)*% and Lj := (—1)'' which take values 
±1, the energy can be rewritten in the form 

^I,r,s = - E ['JR^iS^S^+l + 1) + KL,] . (19) 
j 

In the sequel we will use both the indices r, s, / and R, S, L where no confusion arises. The 
states |^'i,r,s) have 3A^ binary indices and hence the number of such states is exactly equal to 
the dimension of the Hilbert space, hence they will comprise the full energy spectrum, provided 
that we can show they are independent. To investigate this question, let us look at the inner 
product of these states. 

Using (ITSl) . and (IT^ we find 

(^l,r,s|^l',r',s') = ^(^+|-BlAr',s'Ar',s'-Bl/|^7+) 

= 5i^i/(fi+|i3iAr+r',s+s'|^+)- (20) 

We now note that the operators Ar+r',s+s' generate only open or a homologically trivial loop, 
while the operator Bi only generate homologically trivial loops, and hence the above matrix 
element vanishes unless (r, s) = (r',s'). Moreover by expanding Bi it is clearly seen that 
{Q^\Bi\Q^) = 1. Hence we find 

(^l,r,s|^l',r',s') = <^l,l'^r,r'5s,s' 

The independence of these states and the equality of their number with the dimension of 
Hilbert space, indicates that they comprise the complete spectrum of the Hamiltonian. It is also 
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instructive to note the symmetry of the spectrum. If we indicate by s the binary complement 
of the indices s (i.e. Sj = 1 + Sj) and similarly for other indices, we find from (ITS!) that 



Ej-MK) = E,,r,s(-J,-i^). (21) 

The first relation expresses the two-fold degeneracy which is the result of the topology of the 
surface, that is the action of the cycle Wz on any state, produces another state with the same 
energy. The second relation indicates how the spectrum is affected if we invert the coupling 
constants J and K around 0; the spectrum should be inverted around the values Vj = Ij = 0. 
One should also note that only the ground and the top states have two-fold degeneracy and 
the degeneracy of the other states is much larger. 



2.3 The partition function and averages of different observables 

From the complete spectrum it is straightforward to calculate the partition function. One 
writes 

Z{(3, J, K) = tr(e-^^) = ^ e-^^^-^-^ 

= Ene^'''^^''^-^''^-''^Ene^''^^=:^o(/5,J)^i(/?,i^), (22) 

R,S j L j 

where the last equality defines the partition functions Zq := Zq{(3, J) and Zi := Zi{[3, K). It is 
obvious that Zi = 2^ cosh^ (3K. Using a transfer matrix, we find 

^0 = E n 2 cosh (3 J{Sj^iSj + 1) = 22^(cosh2^ f3J + sinh^^ (3J). (23) 
s 

The full partition function is therefore given by 

Z{(3, J, K) = 2=^^(cosh2^ (3J + sinh^^ (3J) cosh^ /3K. (24) 

In the thermodynamic limit, the average energy is obtained from the partition function to be: 
E 8 

(— ) = ——\nZ(f],J,K) = -{2Jtanhf3J+ Ktanh(3K). (25) 

In the same limit the entropy is found from = (1 — P^) In Z to be 

S 

— = 3 In 2 + 2 In cosh /3 J + In cosh PK - 2/3 J tanh 2/? J - tanh ^K. (26) 

The entropy can be written as a sum of two terms, namely S = N{Ss{(3J) + Sp((3K)), where the 
first one is the contribution of the vertex terms and the other is the contribution of plaquette 
terms, and 

Ss(x) = 2 In 2 + 2 In cosh x — 2x tanh 2x 
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Figure 2: Color online. The entropy of the ladder as a function of the couplings and tempera- 
ture. 



Sp{x) = ln2 + Incoshx — xtanhx. (27) 

Figure ([2]) shows the total entropy as a function of /3J and PK. 

An interesting non-local observable is the string operator defined in One finds 

(W.) = z^p\k) E(^i|Ar,siy.Ar,s|vl/i)e-^^'.--. (28) 
Passing Ar,s through W^, one finds 

(vI/i|Ar,,W^,Ar,s|v&l) = (-l)^'i+'^^+--, (29) 

where we have used the commutation relations of the operators and also the fact that {W^Bj = 
BjWx Vj) and consequently W2:|\E'i) = I^E'i). Inserting this into (l28l) one arrives at 

m = ^ E n R,e^'^^^'^-^'^^'\ (30) 
^0 R,s j 



where again an appropriate transfer matrix gives the final result 

1 sinh^ 2/3 J 

2^cosh2^/5J + sinh^^^ (3 J 



(^^) = ^ ....2..., 1.2... (31) 
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In the thermodynamic hmit this gives 

{»-.)={; (32) 

which shows a phase transition at zero temperature. It is easy to see that {Wz) = at all 
temperatures. 

3 Reduced density matrices and entropies of different 
subsystems 

In this section we derive the reduced density matrices of different subsystems at finite tem- 
perature. These subsystems are denoted by A (the spins on one of the legs of the ladder), B 
(the spins of the totality of all the rungs of the ladder), and C (the spins of a subset of the 
rungs), see figure ([3]). The significance of the subsystems B and C is that B corresponds to 
a topologically nontrivial loop in the surface, while C corresponds to a trivial curve and one 
expects that this difference of topology shows itself in the entropy of these subsystems. As we 
will show, this is indeed the case. 

To prepare ourself for the calculation of the reduced density matrix at finite temperatures, in 
each case we first derive the reduced density matrix of the relevant subsystem when the whole 
system is at the state |\E'i). The corresponding density matrices will be denoted by cta, ctb, and 
ac- This will pave the way for determination of the reduced density matrices at arbitrary tem- 
peratures which will be denoted by Pa, Pb and pc- We also use the notation to denote 
the restriction of the product state \ to the subsystem A with similar notations for B and C. 



3.1 Subsystem A: One leg of the ladder 

We have 

■.= tr^„^f-:.j,„mm, (33) 

where for any subset /, / means that we take the trace over the complement of /. To calculate 
the trace, we note that the state |\E'i) can be written as follows 

ivi/,) = -i= y: (-i)S-'^"^^5r---i?nfi+) (34) 

V mi ,m2,---'mjv 

where the indices rrii take the values or 1. We now use the fact that Bi = ZiZi^iZi/Zi« and 
take the trace over the product of the B^ operators to arrive at 

K'^fT-N'^^^'^^" ■ ■ ■ 5n^+)(^+i5f^52^ ■ • • 
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Figure 3: Color online. The links with filled circles comprise the subsystems A, B, C and D. 
At any temperature, the subsystems A, C and D are in maximally mixed states. Only the 
subsystem B which has a non-trivial topology has a different state, Eq. fH9l) . 



where 



From the above equation we find 

i 

which when inserted into (!35|) this gives 

\mi,---mjv 



(35) 
(36) 

(37) 



2N 



(38) 



This shows among other things that at zero temperature, each of the two legs of the ladder 
are in a maximally mixed state. The interesting point is that this situation persists at all 
temperatures. To see this, we take the upper leg as our subsystem A, since in this case, the 
analysis will be greatly simplified. Since the system is symmetric, whatever we obtain will also 
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be valid for the lower leg. We have 

tr^(|vl/,,r,s)(v&i,r,s|) = tr^ ( Ar,s | ^i) (^1 1 Ar,s) = tr^(|^i)(^J) = ^U, (39) 

The reason for taking A to be the upper leg of the ladder is that we could pass through the 
operator Ar,s cyclically within the trace and arrive at the simple result that Pa{T) = i^Ia- 
Obviously any subsystem of A will also be in a maximally mixed state. 

3.2 Subsystem B: All the rungs of the ladder 

Consider now subsystem B, the full set of rungs of the ladder. We first derive the reduced 
density matrix ctb, when the whole ladder is in the state \^\). Using the decomposition ( l34l) 
and the structure of the Bi operators ([2]), we find 

<y, = ^ E (-i)^^''^'"^+''^i^(fi+|iVm,k|f^+)J$„.)(<fk| (40) 

where 

Nrn,^.■■=X{{Z,Ze^^^^^ (41) 

i 

and 

1$^) := {Z^Z^r- ■ ■ ■ {ZNZ^r^\n+)^. (42) 

Using the fact that A^m^k|^+)^ = (^m.k, and inserting the result in (HUI) . noting the two- 

to-one correspondence between the indices rrii and the powers of Zi and rearranging terms, we 
obtain 

^B = ^ E (43) 

gi.-'-'jjv-i 

where 

|<|.q) ■.= ZrZf---Z'^^'''--'^-'\n+)B (44) 

Note that this is independent of the index set 1 of the state I^E'i). Also in each state I'I'q) the 
flip operators come in pair, so this state is an even parity state, i.e. a state where an even 
number of spins have been flipped from + to -. The state cr^ is thus a uniform mixture of all 
even parity states. Call this density matrix a^^^"' , Let us now consider finite temperatures, for 
which we have to calculate 

^^,5;^jv(l^l.r,s)(^l,r,s|) = tr^_5;^^(Ar,s|^l)(^l|Ar,s) = Arff^^^^^^Ar, (45) 

where Ar = Hi . The reduced density matrix at finite temperature will now be given by 

PsiT) = ^ E e-^^^-Ka'^'^'K. (46) 

^0 r,s 
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From the above definition of Ar, and that of a^^'^"' and |$q), one finds that 

; I'riLi^ (47) 

where the uniform mixture of odd-parity states and |r| denotes the degree of r, i.e. 

|r| = n + r2 H Tat. 

Inserting this into (146!) yields 

PBiT) = E ^'^"'-^ ((1 + (-1)1^1)^'=''^" + (1 - (-1)1^1)^°'^'^) . (48) 

^^0 r,s 

Using the expression of the parity |r| and Eqs. (1251) and (12^ we find the following simple 
expression 

PB(r) = ^(l + (iy,))a^^'=" + ^(l-(t^,))a'"^'^ 

= ^/b + ^(W,)(^'^^^"-0, (49) 

where use has been made of the fact that \{a'^'"^'^ + a°'^'^) = ^Ib- The entropy of this state, 
which is a mixture of orthogonal states, can now be readily calculated. A straightforward 
calculation gives 

^(pB) = Ar-l + i/(i±|^), (50) 

where H{p) = — ploggP — [1 — p) log2(l — p), is the Shannon entropy function. 

Figure (jl]) shows Sp^ — {N — 1) for several values of system sizes A^. In the thermodynamic 
limit, there is a sharp rise in this quantity only at zero temperature, but for finite A^, it is also 
seen that there is an almost sharp rise at finite temperatures. 



3.3 Subsystem C: A subset of the rungs 

Once the density matrix of the subsystem B is obtained, we can trace out any number of the 
rungs to find the reduced density matrix of the remaining subset of rungs. Using (H9|) and 
noting that on taking the trace over any subsystem, the contribution of a'^^'^"' and a""^"^ cancel 
each other, one arrives at the simple result that 

Pc = ^\Ic, V T. (51) 

where C is any proper subset of the rungs and |C| is the size of this subset. It is interesting to 
note that how the topology of the surface is reflected in the entropy of its subsystems. 
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Figure 4: Color online. Thermal entanglement of subsystem B with the rest of the lattice as a 
function of temperature for different system sizes. The system sizes are 10, 20, 30, 40 and 100. 

3.4 Thermal entanglement of two spins 

From the reduced density matrices found in previous subsections, we know that any two spins 
on a leg of the ladder are in a maximally mixed state p = |/ and hence there is no thermal 
entanglement between such spins. The same is true between any two spins on the rungs of the 
ladder. In fact it has been shown [27J that in the ground state, there is no entanglement between 
any two qubits. However one can ask if at higher temperature some degree of entanglement 
is caused by thermal fluctuations. This indeed happens in some spin systems, below a certain 
threshold temperature. To investigate this, we compute the reduced density matrices of the two 
spins, say 1' and 1" in figure ([3]) on the two legs, opposite to each other. Call this subsystem 
D. The first step is to calculate tT'_|'0i^r,s)(V^i,r,s| which is equal to 

^tr^ (Ar,s5i|i^+)(n+|Ar,s5i) = ^ZtHr^ iB,\n^){n^\B,) Z^l 

. In calculating the trace, one can use the cyclic property of the trace and move around all the 
terms (1 + {—lY'Bi) except the term 1 + {—lY^Bi (which acts nontrivially on the space D) and 
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use the property (1 + (— = 2(1 + {—lY'Bi) which after some algebra gives 
tr^ {B,\n^){n+\B,) = (^{1 + (-1)''S.)(1 + (-l)'^5i)|fi+)(^^+|(l + (-l)'^5i)j 

= 2^-4r. ((1 + {-iy^B,)\n^){n^\{i + i-iy^B,)) 

= 2^-Ml + +)(+ + | + | --)(-- |)^,^„, (52) 

where in the second hne we have used the fact that the closed loops generated by any of 5j's 
cannot be compensated by Bi to make a non-vanishing trace and the third line is the result of 
explicit expansion and calculation. We will then have 

P. = ^ E e-'^^-^^/d + +)(+ + I + I --)(-- (53) 

l,r,s 

Acting the operators on both sides and performing the above simple calculation with the 
help of the transfer matrix, we find that po = \Id, which means that there is no thermal 
entanglement between these two spins. 



4 The complete spectrum of the three-leg ladder 

What has been done for the two- leg ladder can be extended to three- leg ladder, Fig. ([5]), without 
much effort. The first step is to define the states 

1^1) = :^ 11(1 + (54) 

and then the trick is to find a suitable generalization for the operators Ar.s, so that their action 
on the above state, produces the correct number of independent eigenstates of the vertex 
operators. The suitable generalization is as follows 

N 

Ar,s,t ^=EzrzpZ^ (55) 

i=l 

where the flipping operators correspond to the links shown in Fig. 

One can now easily verify the following commutation relations, where we use Aj', A'^ and Aj 
for vertex operators on site j for the lower, middle and upper legs of the ladder respectively: 

AjAr,s,t = (-l)^^Ar,s,tA+ 
^°Ar,s,t = (-l)^-^+^^+^'^+*^Ar,s,tA; 

AjAr,s,t = {-iy^Ar,s,tA-. (56) 
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Figure 5: The labeling used in the text for the links on the 3- leg ladder. 



A 



B 

Figure 6: The flipping operators are chosen form the set of links shown in bold, A) for the 
two-leg ladder, B) for the three-leg ladder. 

One can proceed along the same way as detailed in section ([2]) for the two-leg ladder and 
show that the states |\l/i,r,s,t) := Ar,s,t|^i) are energy eigenstates with energies given by 

E = -JY,{RiT,Si_iS, + Ri + Ti) (57) 

i i 

where again we have used Ising-like variables, i.e. Tj := (—1)** instead of the binary variables. 
Moreover the number of these states is 2^^ which is equal to the dimension of the Hilbert space 
and they are orthogonal. (See the reasoning following Eq. fl20|) ). 

The basic point is that every conceivable combinations of the flipping operators in Ar,s,t, 
corresponding to Fig. ([6]) when acting on will produce only open or homologically trivial 

loops of negative spins which will certainly be orthogonal to the state generated by B\ acting 
on The partition function turns out to be 



Z(/5, J, K) = 2^^' cosh"'" [cos\f'' f3J + sinh'^' pj\ (58) 

5 Discussion 

We have determined the complete spectrum of the Kitaev model on a spin ladder and from 
there we have determined the reduced density matrices for its various subsystems at finite 
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temperature. We have shown that on two and three-leg ladders, the model is equivalent to 
particular types of one dimensional classical Ising models, models with different spins on the 
sites and links. 



To what extent this study can be pursued for the two dimensional lattice, ( i.e. a torus)? 
On a lattice with A^^ sites, and 2N'^ links, the Hilbert space dimension is 2^^ . We can already 
construct 2^^ (un-normalized) states of the form |\E'i) := Ylpi)- + Bp)\VLj^) which are en- 

ergy eigenstates. To find more states, we have to find subsets / of flipping operators and then 
construct operators of the form Ag := Hie/ ^t^ ^'^d energy eigenstates as |^E's,i) = As-Bi|\E'i). The 
subset / should have the following important property: when acting on no combination 

of links in / should be able to generate a homologically trivial loop of negative spins on the 
lattice. Let Imax be a maximal set of this type with |/max| elements. Then the total number of 
independent energy eigenstates found in this way is 2^ 



For the square lattice of A^^ sites, one such set is shown in Fig. ([7]), where \Imax\ = N"^ ■ In 
fact Imax is nothing but a one-cycle which goes back and forth around the torus, but does not 
wrap around it, and comprises half of the links on the network. We call the maximal cycle, 
since it is the cycle which contains the maximal set of links (the addition of any link to this 
cycle will make a trivial loop out of it), or the excitation curve, since flipping operators chosen 
from it, and acting on |\E'i) create all the excited states. 

Since \Imax\ = N"^ 1 the states constructed as |\E'm,i) = ^mBi\^i) form the whole set of 
energy eigenstates. Numbering the links along the bold-face curve shown in figure ([7]), in 
a consecutive way from 1 to A^^, shows that the energy of such a state is equal to i?s,i = 
— JY.f=i SiSi+i — K Y,i Li, where we have used the Ising type labels Si := (—1)''' and Lj = (—1)'' 
instead of the binary labels Sj and k. This leads to the partition function 



^square lattice(/?, J, K) = 2^^' COsh^' (3K [cosh^^' (3 J + siuh^^' (3J\ (59) 

Remark: We could have taken just such a canonical curve for the two and three leg lad- 
ders, instead of the ones shown in figure ([6]), although it may have rendered the calculations of 
reduced density matrices in these simple cases unnecessarily involved. The curves shown in fig 
^ can be obtained from this canonical curve by the moving up and down the horizontal links 
appropriately which amounts to application of plaquette operators. Such plaquette operators 
affect the eigenstate by only a phase. 

Knowing the full spectrum in this way, will enable us to study the entanglement and many 
other properties of the Kitaev model in detail. Furthermore this knowledge may be useful in 
other more detailed studies of the Kitaev model, where the dynamics of the model is required. 
For example in the study of dynamics of classical and quantum phase transitions in systems 
with topological order, mentioned in the introduction [30]. Another interesting context is the 
study of auto-correlation times of the toric code, which is related to the important problem of 



16 



1 

N 

N-1 



3 

2 

12 3 N-1 N 1 

Figure 7: Color online. The excitation curve for the two dimensional lattice, the curve from 
which the flipping operators are chosen for creation of the energy eigenstates. The links in red 
(the ones in the top and right of the lattice) indicate identical links on the other sides of the 
lattice, due to the torus topology. 

how long quantum information can be protected in topological degrees of freedom in a back- 
ground of inevitable thermal fluctuations. This later problem was first addressed in [321 [33] in 
which among other things, a mapping of the spectrum of the Kitaev model to two uncoupled 
Ising chains were found. 

Another important problem which may be treated in an alternative way by the characteri- 
zation of spectrum in the way shown in this paper, is the problem of sustainment of topological 
order at finite temperatures. This problem has been studied in a number of works. [23] [Mj 135] 
using a different description of the spectrum. 
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